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ABSTRACT 

We investigate the stability of super Earth atmospheres around M stars using a 7-parameter, analytical frame- 
work. We construct stability diagrams in the parameter space of exoplanetary radius versus semi-major axis 
and elucidate the regions in which the atmospheres are stable against the condensation of their major con- 
stituents, out of the gas phase, on their permanent nightside hemispheres. We find that super Earth atmospheres 
which are nitrogen-dominated ("Earth-like") occupy a smaller region of allowed parameter space, compared to 
hydrogen-dominated atmospheres, because of the dual effects of diminished advection and enhanced radiative 
cooling. Furthermore, some super Earths which reside within the habitable zones of M stars may not pos- 
sess stable atmospheres, depending on the mean molecular weight and infrared photospheric pressure of their 
atmospheres. We apply our stability diagrams to GJ 436b and GJ 1214b, and demonstrate that atmospheric 
compositions with high mean molecular weights are disfavoured if these exoplanets possess solid surfaces and 
shallow atmospheres. Finally, we construct stability diagrams tailored to the Kepler dataset, for G and K stars, 
and predict that about half of the exoplanet candidates are expected to habour stable atmospheres if Earth-like 
conditions are assumed. We include 55 Cancri e and CoRoT-7b in our stability diagram for G stars. 

Subject headings: planets and satellites: atmospheres 



1. INTRODUCTION 

The lower effective temperatures of M stars (or red 
dwarfs) — which comprise about three-quarters of the stel- 
lar population — render the detection of orbiting, close-in, 
Earth-like exoplanet s amenable to both the transit and radial 
velocity techniques (Charbonneau 2009). The immense in- 
terest in hunting for exoplanets around these stars stems from 
the fact that the habitable zone is located at ~ 0.01-0.1 AU , 
rather than ~ 1 AU, away from these stars dTarter et al.l2 007). 
A consequence of the close proximity is that the exoplanets 
are expected to be tidally locked (or spin synchronized), 4 
meaning that they possess permanent dayside and nightside 
hemispheres. This expectation has led to theoretical concerns 
that the atmospheres may undergo collapse, because if heat is 
not redistributed efficiently from the dayside to the nightside 
it may lead to the dominant chemical species condensing out 
on the nightside (e.g., Joshi, Haberle & Reynolds 1997; Joshi 
2003). For example, Wordswort h et al.l (1201 ll) constructed 
three-dimensional, general circulation models tailored to the 
study of Gliese 58 Id and elucidated scenarios for atmospheric 
stability. Generally, the discovery of close-in super Earths 
(e.g., Mayor et al. 2009; Charbonneau et al. 2009; Vogt et al. 
2010; Anglada -Escude et al.ll2012t iDelfosse et alj|2012l) has 
inspired a numb er of theoretical studies of their possible at- 
mosp he res (e.g., Merlis & Schneider 2010; Rogers & Seager 
20 lOt ICas tan & M enoul 120111: [Heng, Menou & Phillipps 
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which either focus on specific case studies or explore a 
limited range of parameter space. 
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Strictly speaking, the terms "tidally locked" and "spin synchro- 
nized" a re only synonymous for an exoplanet residing on a circular orbit 
('Heller. Leconte & Barnes 2011). We use these terms interchangeably as we 
do not account for the effects of an eccentric orbit. 



Specifically, we consider an exoplanet possessing a solid 
surface which is enveloped by a shallow atmosphere. Whether 
the surface is made of rock, ice, etc, is unimportant, insofar as 
it provides a solid, lower boundary for the atmosphere. The 
atmosphere is assumed to be shallow enough that it may be 
characterized by advective and radiative time scales which 
are constant over its entire depth. When the atmosphere is 
predominantly radiative, heat is thermally re-emitted before it 
is able to be advected to the permanent nightside. The night- 
side becomes arbitrarily cold, 5 thus leading to the conden- 
sation of its constituents from the gas to the liquid or solid 
phases. When this situation occurs, the atmosphere is unlikely 
to survive for a duration comparable to the stellar age and we 
term the atmosphere to have "collapsed". This situation is dis- 
tinct from the phenomenon of atmospheric escape (e.g., Yelle 
2004; Mur ray-Clay, Chiang & M urray 2009), which concerns 
the (often slow) vertical flow of material out of the gravi- 
tational potential of the exoplanet. When the atmosphere is 
predominantly advective, the efficient redistribution of heat to 
the nightside provides a necessary but insufficient condition 
for the nightside atmosphere to exist in its gaseous form — 
whether it does depends on the thermodynamics of phase 
changes of the chemical species being considered, which pro- 
vides the sufficient condition. If the atmosphere remains in 
the gas phase, then we term it to be "stable". In other words, 
an atmosphere which is sufficiently advective and warm will 
survive for a long time. 

In the present paper, the main question we are addressing 
is: what is the simplest model one can construct to broadly 
understand the stability of super Earth atmospheres? To 
this end, we introduce an analytical framework which allows 
for an efficient, broad exploration of parameter space, while 
keeping the number of free parameters involved to a mini- 
mum. The key outcome of our study is the construction of a 
stability diagram, which allows one to judge if a super Earth 
atmosphere is likely to survive for a duration comparable to 
the age of its host star. In ^2] we describe our methodology. 

5 But see discussion of caveats in i|4. 1 1 
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In Sj3] we present the basic stability diagram as well as varia- 
tions of it applied to GJ 436b and GJ 1214b. In SJU we discuss 
the caveats associated with our model, the implications of our 
results and also apply our stability diagrams to the Kepler 
dataset of exoplanets and exoplanetary candidates, as well as 
55 Cancri e and CoRoT-7b. 

2. METHODOLOGY 

To efficiently explore a wide range of parameter space, we 
construct a 7-parameter model with the following inputs: 

1. The spatial separation between the exoplanet and its 
host star (a); 

2. The radius of the exoplanet {R); 

3. The bulk mass density of the exoplanet {po); 

4. The pressure level associated with the infrared photo- 
sphere (P); 

5. The mean molecular weight of atmospheric 
atoms/molecules (/z); 

6. The number of degrees of freedom of the atmospheric 
gas (naof); 

7. The Bond albedo of the exoplanetary atmosphere (A). 

These basic parameters yield all of the other secondary pa- 
rameters. We assume our exoplanets to be spherical, such that 
the mass is described by M — AnpoR 3 /3. The surface grav- 
ity is given by g = AttGpqR/3, where G is the universal grav- 
itational constant. In addition, the stellar mass (M*), radius 
(i?*) and effective temperature (T*) need to be specified. 

2.1. Stellar Irradiation and Orbital Parameters 
The incident flux impinging upon the substellar point is 

Fo = CsBTta, W 

where ctsb is the Stefan-Boltzmann constant and the irradia- 
tion temperature is 



T;,. 



T* 



R \ 1 / 2 

^) (1-A) 1/4 
a J 



(2) 



If one adopts parameter values appropriate to Earth (and A = 
0), one obtains the solar constant, To « 1370 W m -2 . 

Since the assumption of tidal locking is made, the orbital 
and rotational frequency are equal and given by (assuming 
that M* > M) 

It is thus clear that T m (or Tq) and $7 need to be varied self- 
consistently as one changes a. 

2.2. Thermodynamics 

The mean molecular mass m = /imn (with toh denoting 
the mass of a hydrogen atom) and the number of degrees of 
freedom of the atmospheric gas ndof collectively describe the 
thermodynamics and are manifested macroscopically via the 
adiabatic gas index, 

r = i + — , (4) 

n do f 
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FIG. 1. — Top panel: approximate phase diagrams of various molecular 
species as computed using the Clausius-Clapeyron equation. At a given tem- 
perature, chemical species existing at a pressure above their respective satu- 
ration pressure curves are in the liquid or solid phase. Bottom panel: conden- 
sation temperature as a function of pressure. At a given pressure, chemical 
species existing at a temperature above T CO n are in the gas phase. 
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as well as the adiabatic coefficient ( Pierrehu mbertl20 Tol) . 

_ 2 K 

2 + ndof cp ' 

where 1Z = TV /fi is the specific gas constant and 1Z* = 
8314.5 J K _1 kg -1 is the universal gas constant. The sound 
speed is c s — {Tk^T^/m) 1 / 2 , where fcs is the Boltzmann 
constant. Furthermore, the specific heat capacity, 



c P = 



C2 + n doi )1l* 



2// 



(6) 



is used in evaluating the radiative time scale i ra d. 

As an example, atmospheres dominated by molecular hy- 
drogen have p — 2 and ndof = 5, such that T = 7/5, 
K = 2/7 and 11 = 4157.25 J K _1 kg" 1 . Dry, terrestrial 
air (dominated by nitrogen and oxygen) has a mean molecu- 
lar weight ofpw 28.97, such that 11 « 287 J K" 1 kg" 1 . As 
a furthe r example, we note that the model with n sw 17 from 
Menoy (1201 2t) is termed "Water", with relevance to a possi- 
ble subclass of H2Q-dominated super Earths known as "wa- 
ter wor lds" (Charbo nneau et al.ll2009llRogers & Seagerl2010l 
Bean et al. 201 1). In practice, the effective number of degrees 
of freedom nn n f is no t an integer (e.g., see Table 2.1, pg. 92, 
of Pierrehumbert 2010). For example, CO2 has ndof ~ 6.8, 
O2 has ndof ~ 5.1 and H 2 has ndof ~ 5.2. For simplicity 
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TABLE 1 
Table of Thermodynamic Quantities for Various Molecular Species 



Species 


M 


LavgGJkg- 1 ) 


3con,0 (K) 


T clit (K) 


P sat ,o (bar) 


-Pcrit (bar) 


H 2 


18 


2.7 x 10 6 


5845 


647.1 


1.2 X 10 7 


221 


CH 4 


16 


5.7 x 10 5 


1097 


190.44 


2.1 x 10 4 


45.96 


C0 2 


44 


5.0 x 10 5 


2646 


304.2 


1.1 x 10 6 


73.825 


N 2 


28 


2.3 x 10 5 


775 


126.2 


2.7 x 10 4 


34.0 


2 


32 


2.5 x 10 5 


962 


154.54 


7.4 x 10 4 


50.43 


NH 3 


17 


1.8 x 10 6 


3680 


405.5 


9.2 x 10 6 


112.8 


H 2 


2 


4.5 x 10 5 


108 


33.2 


1.7 x 10 2 


12.98 



and to not obscure the salient physics of the model, we adopt 
n dof = 5. 

To determine if a given chemical species exists in the gas 
or liquid/solid phase, one needs to have knowledge of the sat- 
uration pressure P sa t. If the atmospheric pressure exceeds 
the saturation pressure, then the molecule (or atom) will con- 
dense out of the gas phase and into the liquid/solid phase. 
The dependence of P sat on the temperature T is governed by 
the C lausius-Clapeyron equation (e.g., §2.6 of lPierrehumbertl 
2010), 



dP Rl 



dT 



•*sat-* con,0 
j^2 ' 



(7) 



where T con = L avg /lZ is the characteristic condensation 
temperature. Depending on whether one is describing the gas- 
liquid or gas-solid phase change, one needs to specify the la- 
tent heat of condensation or sublimation, respectively. For 
most gases, the difference between these two latent heats is 
small (w 5-10%); an exception is CO2, where the difference 
is about 33%. For simplicity, we adopt the a verage value of 
these t wo latent heats (L avg ) using Table 2. 1 of Pierrehumbert 
(2010). Furthermore, the value of L avg remains fairly con- 
stant unless T ^> T con .o- Thus, the Clausius-Clapeyron equa- 
tion can be integrated to obtain 



P S! 



P 



sat,o exp 



T, 



con,0 



(8) 



where the normalization constant P sa t,o is computed using the 
pressure and temperature associated with the triple point of a 
given chemical species. In Table [T] we list £ avg , T con and 
Psat,o f or several molecular species. The top panel of Figure[T] 
shows the corresponding (and approximate) phase diagrams, 
which include the triple points associated with each species. 

Our use of the saturation pressure curves is valid as long as 
one is not above the temperature T cr ;t and the pressure P cr it 
associated with the critical point. In this regime, the atmo- 
sphere exists as a supercritical fluid where distinct liquid and 
gas phases do not exist. If the atmosphere is above the criti- 
cal temperature but below the critical pressure, then it exists 
as a gas. Since we typically have P cr it > 10 bar (which is 
larger than the pressures we consider), we assume our model 
atmospheres to be gaseous when T > T cr i t . 

At the infrared photosphere (located at a pressure P), the 
atmosphere remains in gaseous form if P < P sa t> which may 
be recast as the condition T > T cou where 



T, 



con.O 



ln(P sa t,o/P) 



(9) 



is the condensation temperature. The logarithmic factor 
ln(P sat /P) implies that T con has a somewhat gentle depen- 
dence on chemical composition. In the bottom panel of Fig- 
ure Q] we see that T con ~ 100 K. An exception is molecular 



hydrogen for which T con ~ 10 K over a wide range of photo- 
spheric pressures. Generally, it becomes less difficult to con- 
dense out a given chemical species as the pressure increases. 

2.3. Infrared Photospheric Pressure 

We assume a thin atmosphere with a Bond albedo A, such 
that a fraction (1 — A) of the stellar flux reaches the surface, 
which resides at a pressure P . The starlight is absorbed and 
re-emitted in the infrared at a photospheric pressure P < P . 
Furthermore, we assume P ~ Pq such that we may apply a 
shallow water model to represent the entire atmosphere be- 
tween P and Po. The atmosphere is also thin in the sense that 
the radiative time scale t ra( j is assumed to be approximately 
constant throughout, at least near the substellar point. We re- 
gard P as a free parameter which characterizes the thickness 
of the atmosphere. 6 

While we model our atmospheres as consisting purely of 
molecular hydrogen or nitrogen in our stability diagrams, the 
assumption is that they also contain minor amounts of other 
chemical species which absorb strongly in the infrared but 
constitute an insignificant fraction of the total mass of the at- 
mosphere. This enables us to define a "mean gas" with a mean 
molecular weight and to compute quantities such as cp. We 
further assume that these trace constituents are scarce enough 
that their partial pressures are miniscule, such that they always 
remain in the gas phase. In this case, atmospheric stability is 
dictated by the main constituent (i.e., H2 or N2) existing in the 
gaseous regime of the phase diagram. The typical collisional 
time between the major and minor constituents of the atmo- 
sphere is assumed to be the shortest time scale in the system, 
such that the atmosphere may be modeled as a single fluid. 
We neglect any effec ts associated with clouds or hazes (e.g., 
see lHengetaT]|2012h . 

2.4. Hydrodynamic Shallow Water Model 

We utilize the linear, analytical, steady-state, two- 
dimensional version of the shallow water model of 
Showman & Polvani (1201 1|) to model t he photosphere of a 
super Earth. (See also IShowman & Po lvani 2010.) The 
shallow water system has traditionally been used to model 
both the terrestrial atmosphere and ocean (Matsuno 1966; 
Longuet-Higgins 1968; Gill 1980), but has been generalized 
by IShowman & Polvanil d2011l) to include radiative forcing 
and drag within the conte xt of spin-synchronized exoplanets. 
(See also Heng & Spitkovsky 2009 and references therein.) 
Several adjustments are required to mimic the variation of var- 



This situation is distinct from that of hot Jupiters, where one parametrizes 
the optical depths by two broadband opacities (Heng et al. 2012). The radia- 
tive time scale is generally not constant for hot Jupiters and varies with depth. 
The variation of shortwave scattering (i.e., albedo) shifts the photon deposi- 
tion depth to different heights, resulting in varying de grees of hea t redistribu- 
tion between the dayside and nightside hemispheres iHeng 2012). 
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ious atmospheric quantities with the strength of stellar irradi- 
ation. 

Firstly, we need to relate the mean shallow water height H 
to a characteristic height in the atmosphere, which we identify 
as the pressure scale height, 



H 



mg 



(10) 



Secondly, we set the thermal forcing term 5o (which has phys- 
ical units of cm s _1 ) to be H/t , where to = (gH)^ 1 / 4 ^ 1 / 2 
is the characteristic time scale used to non-dimensionalize the 
forcing and f3 = 20,/ R, and also adopt a zonal wave num- 
ber of 0.5. In other words, the dimensionless forcing is set 
to be unity. The third ingredient is the radiative time scale 
dGoodv & Yung|[l989l) . which is the ratio of the thermal en- 
ergy content per unit area (cpT ur fh where m = P/g is the 
column mass in hydrostatic equilibrium) to the flux To, 



t-rad — 



c P P 



(ID 



which enters into the momentum equations via a radiative 
drag term (-v/t Ia d). 

The fourth — and most uncertain — ingredient is the hy- 
drodynamic drag time scale, which we set to be the rota- 
tional/orbital period, 



_ 2tt 

tdrag - -fi. 



(12) 



On Earth, the sourc e of hydrody namic drag originates from 
the boundary layer (lGarrattlll994l) . which is the transition re- 
gion from a "no slip" to a "free slip" condition between the 
atmosphere and the terrestrial surface. In the Held-Suarez 
benchmark test for Earth, the time scale for R ayleigh drag 
is set to be one Earth day (Held & Suarez| [l994l) . Our simple 
assumption for tdrag is plausible and commensurate with the 
fact that the true source of drag in super Earth atmospheres 
remains to be identified. Even the simplest imp lementation 
of a boundary layer scheme for Eart h dTroen & Mahrtlll986[ 
Friers on. Held & Zurita-Gotorl 120061 ) requires the specifica- 
tion of four additional parameters, which is currently unjus- 
tified given the dearth of astronomical information on this is- 
sue. 

2.5. Basic Trends with Strength of Stellar Irradiation 

To assess the efficiency of heat redistribution from the day- 
side to the nightside hemisphere of a tidally-locked super 
Earth, one needs to compare the rad iat ive to the advective 
time s cale dShowman & Guillotll2002l) . iPerna. Heng & Pontl 
(2012) have shown using three-dimensional simulations of at- 
mospheric circulation — albeit in the context of hot Jupiters — 
that a reasonable approach to computing the advective time 
scale is to evaluate 



^adv 



R 



(13) 



where M max = max{u} and u is the zonal wind speed. 

While we do not expe ct the shallow water model of 
Showman & Polvani (201 1) to reproduce all of the trends ob- 
tained from the three-dimensional simulations, we do expect 
ii mal to generally increase with 2| rr when hydrodynamic drag 
is insignificant. The right panel of Figure [2] confirms this ex- 
pectation for a > 0.1 AU, a basic trend which is supported 
by the results obtained from three-dimensional simulations of 
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FIG. 2. — Variation of the maximum zonal wind speed, non-dimensionalized 
(left panel) and in physical units (right panel), with a (or T; rI .), computed 
using the shallow water model of Showman & Polvani (2011) which is 
specialized to tidally-locked exoplanets. For illustration, we adopt M* = 
O.184M , R+ = 0.203-Rq and T* = 3240 K (appropriate to a M3.5 star) 
as well as P = 0.1 bar and .4 = 0. Shown also are the curves for different 
values of the mean molecular weight fi. 

atmospheric circulation (e.g., Perna, Heng & Pontj |2012l) . At 
a < 0.1 AU, the hydrodynamic drag condition we have im- 
posed to mimic the effects of the boundary layer sets in and 
reverses the trend. We have checked that if t^ rag is kept con- 
stant (and at a large value), then u max increases monotonically 
with decreasing a. The left panel of Figure|2] which shows the 
maximum zonal wind speed normalized by the sound speed, 
illustrates the effects of hydrodynamic drag more succinctly. 

The model also predicts that u max decreases when the mean 
molecular weight increases, because a larger value of fi corre- 
sponds to stronger radiative forcing (i.e., shorter i ra( j). Such 
a trend with j j, is consistent with the simulated results of 
Menou (2012), who — by examining models for GJ 1214b 
with fi « 2, 3 and 17 — found that a higher value of /j, leads 
to a higher ratio of dayside to nightside photospheric flux, 
implying that heat redistribution becomes less efficient as \i 
increases. We have checked that the qualitative trends in our 
stability diagrams are similar whether we adopt tdrag = 2ir/il 
or as a constant, arbitrary value, although the quantitative re- 
sults do differ. 

Three-dimensional numerical solutions of the (more gen- 
eral) hydrodynamic primitive equations (e.g., Vallis 2006) 
predict ~ 1-2 km s _1 zonal flows, which are mildly super- 
sonic, to be present in hypothetical super Earth a tmospheres 
tailored to the study of GJ 1214b (M enoull20T2l) . Our esti- 
mates for the maximum zonal wind speed are broadly con- 
sistent with these calculations. Sup ersonic flows are permit- 
ted in shallow water models (e.g., lAntuonol 1201 0). and the 
assumption of vertical hydrostatic equilibrium does not affect 
the s peed of the horizontally propag ating sound waves in gen- 
eral (JMacDonald. Lee & XiefooOOh . . 

2.6. Inefficiency of Thermal Conduction 

We assume that the rocky cores of our model super Earths 
do not efficiently conduct enough heat from the irradiated 
dayside hemisphere to the cold nightside hemisphere to pre- 
vent atmospheric collapse. We will now demonstrate that this 
is a plausible assumption. 

In one dimension, the conduction equation reads 



dT 
~dt 



a 



d 2 T 
dx 2 ' 



(14) 
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where T denotes the temperature, t the time, x the spatial 
coordinate and 

a = (15) 

Pocp 

is the thermal diffusivity of the material of the rocky core. 
Here, cp denotes the specific heat capacity, at constant pres- 
sure, of the rocky core (and not the atmosphere). Assuming 
that T — T(x, t) — T x (x)T t (t) is a separable function of x 
and t yields 

-Ca T t , (16) 



dt 



which has the solution 



T t = T tQ cxp(-Cani), 



(17) 



where T to = T t (t = 0) and — C is the separation constant. 
We identify C = l/R 2 , such that the characteristic thermal 
conduction time scale is 



t-cond — 



Pqc Po R 2 

4 x 10 7 yr 

R 

R^ 



Po 



( 'p<> 



3 gem" 3 lO^rgK^g- 1 

2 / , \ — 1 



(18) 



k 



cond 



10 10 erg K~ l s- 1 



Rocky material is expected to have a thermal conductivity of 
fccond ~1W mK _1 = 10 10 erg K _1 s _1 and a specific heat 
capacity of cp ~ 10 7 erg K _1 g _1 . Demanding that the 
radiative cooling time is less than the thermal conduction time 
yields the condition 



/A- ,\ 1/3 

rp I "-cond \ 

\vsbRJ 



65 K 



10 10 erg KT 1 s" 



1/3 



R_ 

R(f, 



-1/3 



(19) 



The condition in equation ( fT9l implies that as long as the 
rocky surfaces of close-in super Earths are heated to > 100 
K, thermal conductivity will not operate rapidly enough to 
heat the cold nightside hemisphere. The possibility of signifi- 
cant geothermal heating from the core is not considered in our 
models and remains a topic for future investigation. 

3. STABILITY DIAGRAMS 

3.1. Conditions Related to Advection & Radiation 

Since the atmosphere is predominantly forced by stellar ir- 
radiation, the main effect is the competition between advec- 
tion and radiative cooling. A plausible, necessary condition 
for atmospheric stability is i a( j v < £ m d, which yields 



R< 



c P u n 



<P 



osb3 



W 3 it 3/2 (l-„4)- 3/ V/ 2 . 



(20) 



The maximum zonal wind speed, u max , is obtained from our 
shallow water model (see ij2.41 >. A reasonable approximation 
to the condition in equation ( f20b is to set M m ax = c s , which 
yields 



R< 



Tk B 

m 



1/4 



/ 3c P P 

\4:TtGpq(Tsb 

T ^R- 5/8 (l-A) 



1/2 



(21) 



-5/4 /? -5/8 ri _ ^-5/16 a 5/8_ 



Since free gravity waves have a speed ~ c s , the condition in 
equation (TzTTt may also be interpreted as a constraint on heat 
redistribution mediated by free gravity waves. 

We have made the assumption that if convective instabil- 
ity is triggered, it leads to the redistribution of heat primar- 
ily in the vertical and meridional directions (and not in the 
zonal direction), which is corroborated by three-dimensional 
simulations of atmos pheric circulatio n for tida lly-locked, 
Earth-like exoplanets (lJoshill20 03: Merlis & Schneider 2010; 
iHeng & VodlJOl UlHeng. Frierson & Phillippi|20111). 

3.2. Conditions Related to Thermodynamics of Phase 
Changes 

If advection manages to redistribute heat across the entire 
exoplanetary atmosphere, then it may be described by an equi- 
librium temperature T cq = T- m j \[2. Thus, we may identify 
the temperature in equation (|8) as T — T cq . 

Generally, the atmosphere remains in the gas phase if 

a < max {a con , a cl -it} ■ (22) 

The condensation distance is given by 

R* ( r* 



{\-A) 



1/2 



^ \-^con.O 

while the critical distance is given by 
_R* 

&crit — (-. 



In 



R 



sat,0 

~p~ 



, (23) 



(^)V-4) 1/2 . (24) 

\ -t crit / 



Typically, we have a cl -it < a con because T cr i t > T con . 

In other words, an atmosphere consisting of a single chemi- 
cal species (described by T con and P sa tfl) will remain sta- 
ble if the exoplanet is located sufficiently close to its host 
star, such that the atmospheric molecule (or atom) remains 
in the gas phase. If the atmosphere contains several chem- 
ical species, then P refers to the partial pressure of a given 
species. In the present study, we consider only atmospheres 
consisting predominantly of H2 and N2. 

3.3. Conditions Related to Orbital Circularization 
& Spin Synchronization 

Our application of the Showma n & PolvanH (12011b shal- 
low water model breaks down — and the concerns about at- 
mospheric stability are alleviated or even obviated — if the 
exoplanet is not spin-synchronized, but still (approximately) 
applies if it resides on a mildly eccentric orbit. 7 (See 
Laughlin et al. 2009 for an illustration of the effects of stel- 
lar irradiation on the atmosphere of an exoplanet residing on 
a highly eccentric, but spin-s ynchronized, orbit. See also 
iHeller. Leconte & Barnes! 120 111) The time scale associated 
with spin synchronization is (Bodenhei mer. Lin & Mardling 
12001 

W.(»)(K)(±y (25) 
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where u is the rotational frequency of an exoplanet which 
is not initially spin-synchronized. The time scale associ- 
ated w ith the circularization of the orbit is dGoldreich & Soterl 
[19661) 
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We will see later that the £ a d v = *rad and the t c 
do not intersect. 



(26) 



i* lines typically 



Heng & Kopparla 



In both equations ( 1251 ) and ( 1261 1. the tidal quality factor is 
Q = 27r£ , poa k/Ai?, where .Epeak is the peak tidal en- 
ergy stored and AE is th e amount of energy dissipated 
per forcing cycl e (Knopoff 1964; Goldreich & Soter 1966; 
Bodenheimer, Lin & Mardling 2001). Following Knopoff 
( 1964) and Goldreich & Soter ( 1966), we adopt Q = 10-100 
for rocky exoplanets. Spin synchronization generally occurs 
faster than orbital circularization, unless 



a < R 



14 /uj\ 
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(27) 



If w/O = 100, then the numerical coefficient in equation $T7\ 
is about 17. This condition appears unlikely to be satisfied. 
Demanding that t syn < £* and t c - lrc < £* yield a pair of re- 
spective conditions: 
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A caveat is that if the exoplanet ends up in 1 : N tidal res- 
onances where N > 1, then this tends to homogenize the 
zonal temperature differences, thu s further alleviating co n- 
cerns about atmospheric collapse (Wordsworth et al. 201 1 ). 

3.4. Other Conditions 

In order to retain its atmosphere, a super Earth needs to 
possess a radius which is smaller than its Hill sphere i?H = 
a(M/3M*) 1/3 , which implies that 



a > 



9M* 



1/3 



(29) 



The thermal speed (~ c s ) of the atmospheric gas also needs 
to not exceed the escape speed of the exoplanet, v csc = 
(2GM/i?) 1 / 2 , 
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Equation d30l ) may be regarded as bein g lib- 
eral, as atm ospheric esca pe (e.g., lYellel l2004t 

iMurrav-Clav. Ch iang & Murray] |2009) will still occur 
when v csc > c s . We do not expect the lower limit in equation 
( f30l > to be revised upward by more than an order of magnitude 
if atmospheric escape is properly accounted for (and if it 
occurs at a sufficiently high rate to be relevant over the stellar 
age). 

The habitable zone, which is the range of distances from 
a star where an Earth-like exoplanet can permanently har- 
bor liquid water, may be described by fitting formulae 
based on one-d imensional radiative-convective calculations 
(ISelsis et alJl2007l) : 
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(31) 



where T* = T* 
3.839 x 10 33 erg s 
a 2 = 2.7619 x 10 



5700 K, C* = 4irRla SB T? and £ = 

1 The fitting coefficients are a% = 0.84, 

5 K" 1 , a 3 = 3.8095 x 10~ 9 K" 2 , 
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FIG. 3. — Basic stability diagram in the parameter space of exoplanetary 
radius (R) versus the spatial separation between the exoplanet and its host 
star (a). The hatched regions mark the regimes in which the atmosphere is 
unstable. The pressure associated with the infrared photosphere is set to be 
P = 0.1 bar, while the Bond albedo is .4 = 0. The boundaries of the 
habitable zone are marked by the pair of vertical, triple-dot-dash lines. For 
H2, we have a con fa 24 AU, and thus it is not shown in the diagram. 

a 4 = 1.67, a 5 = 1.3786 x 10~ 4 K" 1 and a 6 = 1.4286 x 10~ 9 
K~ 2 . The quantities ai n and a out represent the inner and outer 
boundaries, respectively. For simplicity, we have adopted a 
cloud-free description of the habitable zone; in such a sce- 
nario, the habitable zone for the Solar System extends from 
about 0.84-1.67 AU. 

3.5. Basic Stability Diagram 

Figure [3] shows the basic stability diagram for super Earth 
atmospheres in the parameter space of R versus a. For the 
condition i syn < £*, we have used w/fi = 100 merely as 
an illustration. As expected, the condition R/c s < £ ra( j is 
less constraining than £ ac i v < ^rad across all of the considered 
values of R, a and [i\ reassuringly, both conditions produce 
similar qualitative trends. The condition Rn > R is over- 
whelmed by both the i a dv < irad and v csc > c s conditions. A 
key feature of Figure [3] is that atmospheres with lower mean 
molecular weights (/1) tend to occupy a larger allowed region 
of parameter space. 

We include in Figure [3] the measured values of R and a 
for GJ 436b and GJ 1214b. For GJ 436b, we use the val- 
ues of R ~ 3 .96.Bm and a w 0.0287 AU as reported by 
Knuts on et alj ( 1201 ll) For GJ 1214b, we use t he values of 
R « 2 .66i? m and a w 0.014 AU as reported bv lCarteret al.l 
d2011l) . Since both GJ 436b and GJ 1214b have a bulk den- 
sity of /?o ~ 2 g cm~ 3 , we adopt this value. We adopt stel- 
lar parameters appropriate to a M3.5 star (M* = 0.184Af Q , 
R+ = 0.203-R©, T* = 3240 K) as a compromise between 
the M2.5 a nd M4.5 host stars of GJ 43 6b and GJ 1214b, re- 
spectively (Ehrenreich & Desert] |2011|) ; we will refine these 
parameter values later in £13.61 We also use t+ = 6 Gyr to be 
consistent with the reported ages of the host stars. 

Another key feature of Figure [3] is that some super Earths 
(R > i?e) residing in the habitable zone (equation OTI ). with 
nitrogen-dominated ("Earth-like") atmospheres (/1 = 28), 
do not possess stable atmospheres. Certainly, this state- 
ment depends on the assumed value of the infrared photo- 
spheric pressure P. Although we do not explicitly include 
it in our stability diagram, we have repeated the exercise 
for oxygen-dominated atmospheres (/1 = 32) and find that 
the conclusions associated with them are similar to those for 
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nitrogen-dominated atmospheres. Thus, it is reasonable to 
te rm nitrogen-dominated atmo spheres as being "Earth-like". 
Kite, Gaidos & Manga (201 1) have identified a pair of cli- 
mate instabilities, related to weathering, which are suppressed 
for atmospheres with high radiative efficiencies (i.e., short 
irad)- Specifically, the enhanced substellar weathering in- 
stability (ESWI) operates when advection is efficient enough 
to counteract the increase in temperature near the substellar 
point, which decreases the rate of weathering. The supply of 
greenhouse gases (via outgassing) now outweighs its removal, 
implying that a runaway greenhouse effect may be triggered. 
Thus, ESWI is activated when advection is efficient, but it 
acts on geological, rather than dynamical, time scales. This 
has the implication that even if irradiation is strong enough to 
suppress advection and the ESWI, atmospheric collapse will 
sti ll occur. At face value, our re sults are orthogonal to those 
of iKite. Ga idos & Manga (2011), implying that super Earths 
which manage to overcome the ESWI need to also possess 
sufficiently advective atmospheres to avoid collapse. Future 
work will elucidate the restricted parameter space in which 
both effects are precluded. 

3.6. Stability Diagrams Applied to GJ 436b and GJ 1214b 

To construct stability diagrams in the parameter space of /x 
versus P, we rewrite the condition in equation I 



H< 



(2 + n doi )lZ* 
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T- 3 R~ 3/2 (1 - A) 
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For GJ 436, we use M„ = 0.45M P , R+ = 0A6R Q and 
T* = 3585 K (lEhrenreich & Desert) 120111 For GJ 1214, 
we use M* = 0.15Mm, R * = O.21i? and T* = 3026 K 
(lEhrenreich & Desertj201 1|) . The surface gravities of GJ 436b 
and GJ 1214b are approximately 13.0 and 9.0 m s~ 2 , respec- 
tively. 

Figure|4]shows the stability diagrams as functions of p and 
P, specialized to GJ 436b and GJ 1214b. In general, atmo- 
spheres with high mean molecular weights (/i) and low P 
values are unlikely to be stable due to the violation of the 
iadv < trad condition. 

The Spitzer transmission spectrum of GJ 436b reveals 
variations indicative of spectral absorption features asso- 
ciated with carbon monoxide and methane (iKnutson et alj 
1201 lb . Our stability diagram for GJ 436b allows for a 
stable atmosphere with /j, ~ 10 if the photosphere re- 
sides at P > 1 bar and A = 0. Non-zero values 
of the Bond albedo allow for higher values of /x. In- 
ferences about the atm ospheric compos i tion o f GJ 1214b 
remain controversial (Rogers & Seag erj | 2010t lBeanetal. 1 
20111: Crolletal. 2011; Des ert et alj 1201 It Ide M ooii et al. 
2012; Miller-Ricci Kempton, Zahnle & Fortnevl 1201 2h . in- 



cluding the possibility that a reduced value of the radius — 
which may obtain from an updated estimate of the distance to 
GJ 1214 via improved parallex measurements — is consistent 
with an atmosphere-less exoplanet. Detecting or ruling out 
the presence of a Rayleigh scattering slope in the spectrum at 
sh orter wavelengths is a valuable atmospheric diagnostic for 
^ dde Mooii et al.ll2012l) . 

4. DISCUSSION 

The theory of (exo)planet formation remains fraught with 
uncertainties. Therefore, an attempt to construct a formation 
model for super Earth atmospheres — which are probably sec- 
ondary (i.e., due to outgassing) and unlikely to reflect the el- 
emental compositions of the primordial nebulae from which 
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FIG. 4. — Stability diagrams in the parameter space of mean molecular 
weight (p) versus infrared photospheric pressure (P), specialized to GJ 436b 
(top panel) and GJ 1214b (bottom panel). The hatched regions mark the 
regimes in which the atmosphere is unstable for a given value of the Bond 
albedo (.4). The stellar (Af*, R* and T*) and exoplanetary (R, a) parameters 
are all set to their observed values for each exoplanet (see text). 
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FIG. 5. — Maximum exoplanetary radius i? max allowed by the t at j v < t ra d 
condition across the spectral types M0 to M9. The left panel shows the set of 
Rmax curves, versus a, for fi = 28 and P = 0.1 bar. The crosses mark the 
location of the center of the habitable zone for each spectral type. The right 
panel shows the i? max value corresponding to the center of each habitable 
zone as a function of the spectral type and for different combinations of the 
mean molecular weight fi and infrared photospheric pressure P. 
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FIG. 6. — Basic stability diagrams applied to the Kepler dataset for G (top 
panel; 1135 objects) and K (bottom panel; 613 objects) stars. The photo- 
spheric pressure is assumed to be P = 1 bar in these diagrams (see text). 
The triple-dot-dash lines denote the boundaries of the habitable zones. For 
the stability diagram applied to G stars, we include 55 Cancri e and CoRoT- 
7b as well as highlighting Kepler-lOb and Kepler-lOc. For the top panel, we 
have a con m 10 and 136 AU for N2 and H2, respectively; for the bottom 
panel, these values are about 5 and 7 1 AU. 

they formed — remains an uncon strained exercise in some in- 
stances (e.g., Migu el et alj|201 ll) . From this perspective, our 
analytical framework may be regarded as a survival — rather 
than formation — mode l, similar to what has bee n constructed 
for planetesimal disks dHeng & Tremainel uOlO). 

4.1. Caveats, Future work & Comparisons to Previous Work 

We have assumed that when advection is inefficient, the 
nightside of a tidally-locked super Earth becomes arbitrar- 
ily cold. In reality, a non-negligible amount of heat may 
be advected to the nightside when £ adv ~ ^rad- The tem- 

S ierature difference betwee n the dayside and the nightside is 
Showman & Guillotll2002l) . 

AT ~ T irr [1 - CX P (-Wtrad)] • < 33 > 

Since AT ~ T nl — r n i g ht, where T n i g ht is the average tem- 
perature of the nightside hemisphere, we obtain 
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such that i adv /i rad -> 00 when T night -> 0; when T night -> 
T irr , we get t adv /t rad -> 0. Adopting T* = 3240 K, R* = 
0.203i? Q , a = 0.001 AU and A = yields t adv ~ 6i rad if 



other words, chemical species with low condensation temper- 
atures (such as H2 with T con ~ 10 K) may remain in gaseous 
form even when the advective exceeds the radiative time scale 
by a factor of a few. Specifically, it is possible for a trickle of 
heat to maintain an atmosphere if it is hydrogen-dominated 
even in cases where i adv < i rad is unfulfilled. Thus, the 
tadv < t rad curve on the stability diagram should be viewed 
not as a sharp boundary, but a smooth transition. In general, 
chemical species with very low condensation temperatures 
probably do not require t adv < t lad as a necessary condition 
for atmospheric stability. The gaseous state of these chemi- 
cal species, with low condensation temperatures, may also be 
maintained by internal heat, which we do not consider. 

Several other caveats provide motivation for future work. 
We have only considered atmospheres consisting of a single 
chemical species, but this approximation should not be too 
severe in situations where the given gas (e.g., H2, N2) is the 
dominant constituent. More refined calculations may consider 
cocktails of greenhouse gases. Another assumption we have 
made is that our model, rocky exoplanets are not covered by 
oceans. Oceanic and atmospheric dynamics couple in a non- 
trivial manner, over a broad range of time scales, with conse- 
quences for the heat redistribution between the day- and night- 
side hemispheres, and are beyond the scope of the present 
study. Our hope is that our simple framework will motivate 
future work on a broad range of atmospheric constituents and 
not just Earth- or even Solar System-centric ones. More re- 
fined trains of thought include the consideration of feedback 
mechanisms such as the water ice-albedo feedback, which has 
been shown to operate weakly on tidally-locked exoplanets 
with Earth-like stellar irradiation and atmosph eric conditions 
dJoshi. Haberle & Revnoldil^lloshillxm 

Finally, we compare the present study to previous studies 
utilizing three-dimensional simulations of atmospheric circu- 
lation. For example, the rocky, CQ2-do minated super Earth 
model considered by IWordsworth et al.l (12011b has i rad ~ 
10 7 -10 8 s, but they do not report their computed zonal wind 
speeds and thus we are unable to determine t adv . However, 
we can infer from their Figure 1 that even though their model, 
in 1:1 tidal resonance, appears to have t adv > t rad , the night- 
side temperatures are about 280-290 K. These temperatures 
lie above the condensation temperature of CO2 at 20 bar, 
which is about 240 K, and thus atmosphe ric collapse does not 
occur. Joshi, Haberle & Reynolds (1997) do not provide the 
value of cp used; if we use cp — 10 7 erg K _1 g _1 , then we 
obtain i rad ~ 10 7 s for their runs with P ~ 1 bar. Since 
"max ~ 10 m s _1 , we have t adv ~ 10 5 -10 6 s (< t rad ) for 
their runs with R = 6400 and 13000 km, implying that these 
model atmospheres are m ostly adv ective, which is co nsistent 
with the findings of Uoshi. Haberle & Reynolds! dl997l) . 

4.2. Maximum Stability Radius 

Our basic stability diagram (Figure |3) demonstrates that 
there is a maximum exoplanetary radius i? max above which 
the (thin) atmospheres of super Earths are unstable. In the left 
panel of Figure|5] we show i? max as a function of a across the 
spectral types M0 to M9. The less intensive irradiation from 
the cooler M stars leads to more lethargic radiative cooling in 
the atmospheres of their super Earths, which increases i? max 
across all values of a. However, this effect is offset by the cen- 
ter of the habitable zone moving inwards, towards the star, for 
later spectral types. The collective result is that i? max is insen- 
sitive to the spectral type. Rather, it is somewhat sensitive to 
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the mean molecular weight (p) and the infrared photospheric 
pressure (P), as shown in the right panel of Figure [5] Atmo- 
spheres with lower photospheric pressures and higher mean 
molecular weights have lower values of i? max , implying that 
the condition for stability is more stringent. For Earth-like at- 
mospheres (p = 28) and g ~ 10 m s~ 2 , the photospheres of 
super Earths with R ~ 1-10P© need to be located at P > 0.1 
bar for stability. 

4.3. Application to Kepler Exoplanets and Exoplanet 
Candidates 

As a final demonstration of the power of our stability dia- 
gram, we apply it t o the Kepler dataset (Borucki et al.ll2.01 it 
Batal haet alJl2012l) . We include all exoplanets and exoplanet 
candidates with R < 6P®. 8 We construct two diagrams for 
G (5200 K < T* < 6000 K) and K (3700 K < T* < 5200 
K) stars using the average values of the stellar parameters for 
each sub-sample. We adopt t* = 5 Gyr for illustration, not- 
ing that this only affects the conditions on circularization and 
spin synchronization. In the absence of empirical constraints, 
we adopt Solar System-centric values for various parameters 
(po = 3 g cm~ 3 , A = 0.3), but note that our estimates are in- 
sensitive to these assumptions. A key, unknown parameter is 
the infrared photospheric pressure. In Figure |6] we set P = 1 
bar. With P = 1 bar, almost all of the 1 135 Kepler exoplan- 
ets and exoplanetary candidates orbiting G stars have stable 
atmospheres if the atmospheres are dominated by molecular 
hydrogen (p = 2), whereas about 44% of them have stable 
atmospheres if they are Earth-like in composition (p = 28). 
For P = 0.1 bar, these percentages become about 88% and 
1 %, respectively. 

We include in the stability diagram for G stars the measured 
R a nd a values for 55 Cancri e ( R » 2.17P m , a « 0.01544 
AU; lDemoryetal1l20TTt iGillion et al. 2012) and CoRoT-7b 
(R « 1.68P®, a « 0.0172 AU: iLegeret al.ll2009l We 
also highlight the exoplanets Kepler- 10b (R « 1.416.R®, 
a « 0.01684 AU; Batal haet al.ll2011l) and Kepler- 10c (R « 



2.227P®, a « 0.2407 AU; iFressin et al]|2011l) . It remains 
possible that 55 Cancri e, CoRoT-7b and Kepler- 10b may 
maintain minimal atmospheres — which do not extend beyond 
the day-night terminators and are dynamically dominated by 
vertical, rather than horizontal, flows — established through 
vapour saturati on equilibrium with their continuously erode d 
rocky surfaces (ICastan & Menoull201 It IValencia et al.ll201 ll) . 

For the K stars (613 exoplanets and exoplanet candidates), 
we obtain similar numbers: for P = 0.1 bar, we get 97% 
(p = 2) and 4% (p = 28); for P = 1 bar, we get 100% (p = 
2) and 70% (p = 28). Curiously, the exoplanet candidates 
and exoplanets residing within the habitable zones of both the 
Kepler G and K stars are likely to possess stable atmospheres, 
unlike the situation with M stars. 

If these exoplanetary atmospheres are Earth-like (P = 1 
bar, p = 28), then about half of them are expected to be sta- 
ble. The continued expansion of the Kepler dataset and the 
eventual, inexorable pursuit of follow-up observations to fur- 
ther characterize these exoplanets will allow our predictions 
to be tested. For example, the photospheric infrared emission 
from atmosphere-less super Earths is likely to be emanated 
directly from their rocky surfaces, such that the dayside to 
nightside flux contrast will appear marked; these exoplanets 
will also exhibit flat transmission spectra. The scrutiny of 
atmosphere-less super Earths may provide a survey of their 
surface compositions (Hu, Ehlmann & Seager 2012). 
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We remain agnostic about the theoretical demarcation between Earth- 
and Neptune-like exoplanets. 
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